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Abstract 

We study the dynamical properties of a two-dimensional ensemble of self- 
propelled dumbbells with only repulsive interactions. After summarizing the 
behavior of the translational and rotational mean-square displacements in 
the homogeneous phase that we established in a previous study, we analyze 
their fluctuations. We study the dependence of the probability distribution 
functions in terms of the Peclet number, describing the relative role of active 
forces and thermal fluctuations, and of particle density. 
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1. Introduction 

Active matter is characterised by the continuous partial conversion of 
internal energy into work. Some examples, at different scales, are the cy- 
toskeleton, bacterial colonies, algae suspensions, bird flocks and schools of 
hsh. Self-propelled units can also be artihcially realized in the laboratc^ in 
different ways, for example, by surface treatment of colloidal particles [l|, l2|. 
All these systems live, or function, in conditions far from thermodynamic 
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equilibrium and pose challenging questions to non-equilibrium statistical me¬ 
chanics. Active matter exhibits non-trivial collective properties that have no 
analogue in passive materials such as large scale coherent motion in the ab¬ 
sence of any attractive interaction and a phase separation into an aggregate 
and a gas-like phase. Several review articles are devoted to this rapidly de¬ 
veloping field of research 0, H B B B B B 0, Q Q Ei- 


The diffusive properties in these systems are of particular interest. A 
number of experimental and numerical studies addressed how the diffusive 
properties are affected by self-propulsion and the density of the suspension; 
some focused on the dynamics of passive tracers immersed in the active 
bath [l^, others focused instead on the mean-square displacement of the 


active particles themselves [1^, [ly, Il7|, [18 
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An interesting model of active matter is one in which the active compo¬ 
nents have the elongated shape of many natural swimmers. A first study 
of the phase diagram of such a system with active dumbbells molecules 
appeared in j^, [^. The analysis of effective temperature ideas j^, 


I, and 

the averaged rotational and translational mean-square displacements were 
presented in and respectively, for a two-dimensional system. 

In this paper we recall some of the results in these publications and we 
extend the analysis to the fluctuations of translational and rotational degrees 
of freedom. In Section 2 the dumbbell model is very briefly explained. In Sec¬ 
tion 3 the numerical results for the translational and rotational fluctuations 
in the interacting active system are presented. A discussion will complete 
the paper in Section 4. 


2. The model 

We briefly present the model and the parameters used in the simulations. 
More details can be found in 23, 24 . The dumbbells are diatomic molecules 


formed by two spherical colloids, elastically linked together via the hnite 
extensible non-linear elastic force 

kr 


Ffene .. / 9 / 2\ ’ 

1 - (rVrg) 

with k > and r = ri — r 2 the vector linking the centres of the spherical 
colloids, with diameter and mass m^. An additional Weeks-Chandler- 
Anderson potential. 


14ca(?^) = 


FLj(r) - I4 j W r < 

0 r > Tr 


2 



















with 


( 2 ) 


VLj(r) = 46 




where e is an energy scale and Vc is the minimum of the Lennard-Jones 
potential, Vc = ‘2}!^ a a-, is added to ensure that the colloids in the same 
molecule do not overlap. The active forces are polar and act along the main 
molecular axis n, are constant in modulus but follow the molecules’ rotations, 
and are the same for the two spheres belonging to the same molecule. 


Fact = Fact fi . (3) 

Fact is directed from the ith colloid (tail) to the i + 1th colloid (head). The 
active forces are applied to all molecules in the sample during all their dy¬ 
namic evolution. We take the interaction between the spheres in different 
dumbbells to be purely repulsive. 

Putting these ingredients together, the dynamic equations are 

2N 

mAii{t) = -'-iYi{t) + Ffene(ri,i+i) + Vi- 

j=0 

mAri+i{t) = --iVi+i(t) - Ffene(ri,i+l) + Vi+l -V ^ + Fact* 

j=0 ^i+l,j 

j^i+l 

with i = 1,3, ...2N - 1, = r* - r^, Vij = |ry| and = 14ca(?^p) with 

Vwca defined in Eq. (|2]). 

The coupling to the thermal bath at temperature T is modeled in the 
manner of Langevin, with 7 the friction coefficient and 77 a Gaussian random 
noise with zero mean, {riia{t)) = 0 , and 

ivMVjbit')) = 2'ykBT6ij6abS{t-t') , ( 5 ) 

with ks the Boltzmann constant, a and b label the coordinates in d di¬ 
mensional space. An effective rotational motion is generated by the random 
torque due to the white noise acting independently on the two beads. 

The surface fraction is 

(j) = N ^ (6) 


^ '^wca p 

dvij r,y 


act 2 ? 


( 4 ) 
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with S'd = T^cr\/2 the area occupied by an individual dumbbell in d = 2, S' 
the total area of the box and N their total number. The spring is supposed 
to be massless and void of surface. We impose periodic boundary conditions 
on the two directions. 

The Peclet number, Pe, is a dimensionless ratio between the advective 
transport rate and the diffusive transport rate. For particle flow one defines 
it as Pe = Lv/D, with L a typical length, v a typical velocity, and D a 
typical diffusion constant. We choose L —)■ ad, n —t Fact/7 and D -)■ = 

kBT/{2'j) of the passive dumbbell to be derived below; then, 

2(TdFact /_x 

ksT 

The active Reynolds number React = ''^d-^act/(c’'d7^) is defined in analogy 
with the hydrodynamic Reynolds number. 


3. Single dumbbell dynamics 


The averaged single dumbbell motion can be derived analytically un¬ 
der the hypothesis that r k, a Details on the calculations can be found 


in [23|, 12^. Here, we simply summarise the main results. Within this ap¬ 


proximation, at absolute times and time-differences that are longer than the 
inertial time-scale tj = md/ 7 , not taking into account its periodic character, 
the angle 9 between the dumbbell’s main molecular axis and an axis fixed 
to the laboratory is a Gaussian random variable with mean (9) = 9q that 
diffuses according to 

( 02 ) = el + 2DRt ( 8 ) 

with 00 the initial angle, t the time-delay, and the angular diffusion constant 


Dr = 


2kBT 

1^1 


(9) 


Averaging over the initial angles, taken from a flat distribution around 0o = 0, 
yields (0) = 0 and, in the long times limit, (0^) —)■ 20^1. In the absence 
of interactions, the angular displacements, A0 between two times that are 
longer than t/ is also Gaussian distributed. If one imposes the periodicity of 
the angles in the interval [—7r,7r] the angular distribution becomes flat. 
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The translational mean-square displacement is ballistic in the limit t 
tj, and crosses over to a very rich behaviour beyond this time-scale, 


(ArD(t) = 4T>cm ^ + 


Fact\" 2 


where 
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1 J Dr 
ksT 


t — 


1 - 


D 
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( 10 ) 


( 11 ) 


is the diffusion constant in the passive limit, F^ct = 0, see [23|. This equation 


presents several time scales and limits. For t/ <C t <C ta = Dn^ one finds 


(ArL) = 4r>J" (+ 


act 


7 


D , 


( 12 ) 


that can still be split into the passive diffusive limit = 4Z1 p^ t for 

tj t < t* ta, and a ballistic regime (Ar^^^) = t* < t ta, 

where the time scale t* is given by 


t* = 



2k bT^ 

act 




(13) 


Note that these two intermediate regimes might be hidden if the system 
parameters are such that t* < tj or t* > ta- They can also be easily confused 
with a super-diffusion behavior (Ar^j^) ~ with 1 < a < 2 if they are not 
well separated [tj ~ ~ to)- In the large Pe limit one has t* <C ta. In the 

last time-lag regime t ta, we recover normal diffusion, 

(ArL) = 4Dx O (14) 


with the diffusion coefficient 


DA{Fact,T,(j) = 0 ) 


kpT ^ , 1 Fact^d V 
27 ^ 2 \ keT J 



(15) 


3.1. The distributions for the single dumbbell 

The infinitesimal increment of the centre of mass position, dxcm, is a 
random variable and its distribution is due to the Gaussian random character 
of the noise and the, for the moment unknown, random character of = 
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cos 9. In order to establish the pdf of the latter, we assume that 6 is uniformly 
distributed in the interval [—tt, vr], pe{9) = (27r)“^. Therefore, for 0 < 6* < tt 





cos 6*) 


1 

2iryi -nj 


(16) 


This calculation can be repeated for —vr < 6 < 0 with the same result and 
the overall numerical pre-factor is hxed by normalisation; 


Pux iP-x) 



for e [-1,1] . 


(17) 


In order to estimate the pdf of the translational displacement of the centre 
of mass position, we now use the over-damped equation for a small increment 
Axcm, with the notation used in 


with Uact 


AXcm = AW + UactAt cos 9 

27 


-^act /7 and AW a gaussian noise, distributed as 


Paw(C) 


1 

- p 2At 


with ( G [—cxo, cxo]. 

The pdf of Axcm is then 


(18) 


(19) 


p( AXcm) 


dCdrZa; 6 



y/AksT'j 

27 


C - fact At n 


Pn^MPAwiO 


The naive comparison of the order of magnitude of the last two terms inside 
the delta-function using ^ ~ (At)^/^ and — 1 yields, after multiplying by 
ad'y/iksTAt), 


Crh \ _ CTdFact 

ksTAt) ksT 


( 20 ) 


In this way the second member corresponds to Pe. Two limits are clear: 

- For Pe left-hand-side in Eq. fl2U]) . so that the active-force induced last 
term in the argument of the Dirac delta can be neglected, the two integrals 
decouple and Axcm is naturally Gaussian distributed. 
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- For Pe S> left-hand-side in Eq. fl20p . so that the noise term in the argument 
of the Dirac delta can be neglected, the two integrals decouple again and 


p(AXcm) OC 


1 



( 21 ) 


with two peaks at Axcm = 't^actAf. 

Otherwise, the double integral yields a complex result: 


p{AXcm) OC 


dn ^exp 


'-1 


( 27 )^ 

2At{4:kBT^) 





( 22 ) 


The displacement between two times. Ax = Xcm(t + to) — Xcmito), is the 
sum of the small increments with each of these independently distributed 
according to Eq. ([22]). 


4. Translational motion at finite density 


The interacting case cannot be solved analytically. We focus here on 
the numerical determination of the centre of mass and angular displacement 
statistics. Details on the numerical method used for solving the dynamical 


equations 


are given m 
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We set rrid = Ud = = e = 1, ro = 1.5, 

A; = 30 and 7 = 10 in proper physical units. These choices assure over¬ 
damped motion and negligible dumbbell vibrations. We used between 15000 
and 20000 dumbbells in each simulation. We fix the strength of the active 
force to be Fact = 0.1 and we vary the temperature in order to access different 
Peclet numbers that we choose to be Pe = 40, 20, 4, 2. The remaining 
parameter is density and we typically use (j) = 0.01, 0.1, 0.3, 0.5, 0.7, see 
Fig.[T]that shows two instantaneous snapshots of the dumbbell configurations 
of systems with 0 = 0.5 (left) and 0 = 0.7 (right). The single-dumbbell 
characteristic time-scales for these parameters are summarised in Table [TJ 


Aspects of the phase diagram and the dynamics of this system were al¬ 
ready established in j^, 23|. At sufficiently low temperature and large ac¬ 

tive force the system phase separates into gas-like spatial regions and clusters 
of agglomerated dumbbells. The dynamic phase transition between homoge¬ 
neous and aggregated phases was determined by the change in behaviour of 
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T 

Pe 

ta 

t* 

0.005 

40 

1000 

10 

0.01 

20 

500 

40 

0.05 

4 

100 

1000 

0.1 

2 

50 

4000 


Table 1: Parameters and characteristic times as defined in Sect.[3l 


the probability distribution function, p, of the local density, 03, [3, 23 

At the critical Pe at which the system starts aggregating the density distri¬ 
bution p not only becomes asymmetric but starts developing a second peak 
at the density of the clusters. Snapshots of typical configurations and their 
analysis along these lines can be found in [2^. In the rest of the paper we 
use sufficiently low Peclet numbers so that the system is in the homogenous 
phase though with important fluctuations, as we will see. 



Figure 1: Typical snapshots of a system with 0 = 0.5, 0.7, linear size 70 Cd and Pe = 40. 
In each dumbbell the green colloid is the head and the red one is the tail, and they are 
joined by a line. 


4-1. Center of mass translational mean-square displacement 

In Fig. [2] we show the center of mass Mean Square Displacement (MSD) 
normalised by time-delay in such a way that a plateau signals normal dif- 












fusion. The two panels display data at T = 0.005 and T = 0.1, under the 
same active force = 0.1. Each panel shows data for eight densities given 
in the key and the case 0 = 0 corresponds to the single dumbbell problem. 
The characteristic times tj, t*, ta (see Table [T]) are shown with small vertical 
arrows. These plots show several interesting features that reproduce, to a 
certain extent, the single particle motion summarised above: 

- In all cases there is a hrst ballistic regime (the dashed segment close to 
the data is a guide-to-the-eye) with a pre-factor that is independent of 0 and 
increases with temperature. (The case t <^tj of the single dumbbell.) 



Figure 2: The center of mass MSD normalised by time-delay, for an active system under 
Fact = 0.1 at r = 0.005, 0.1 (Pe = 40, 2), with different densities given in the key. The 
Peclet number induces a strong qualitative change in (Ar^,^), see the text for a detailed 
discussion. The two dashes in the first panel represent the ballistic dependence ~ The 
dashed segment in the last panel is a guide-to-the-eye for the density dependence of the 
last crossover time-delay that increases weakly with (j). The vertical black arrows indicate 
the single dumbbell time-scales ti and ta-, while the red arrows indicate the single dumbbell 
characteristic time t*, for each case. 


- The dynamics slow down next and, depending on Pe and 0, the normalised 
MSD attains a plateau associated to normal diffusion or decreases, suggest¬ 
ing sub-diffusion. (The case tj t t* ta of the single dumbbell.) For 
instance, there is sub diffusion for Pe = 2 and 0 > 0.4 and Pe = 40 and 
0 > 0.7. 

- Subsequently, the dynamics accelerate with a second super-diffusive regime 
in which the curves for all 0 look approximately parallel and very close to 
ballistic at sufficient high Pe. (The case tj ^ t* ^ t ta of the single 
dumbbell.) 


9 







- Finally, the late normal diffusive regime is reached with all curves saturat¬ 
ing at Da- (The case tj ^ t*, ta ^ t of the single dumbbell.) This regime 
is beyond the time-window explored for 0 = 0.7. 

It is hard to ensure whether the intermediate regime is super-diffusive or 
simply ballistic as the time-scales t* and ta are not sufficiently well separated 
at high Pe and not even ordered as t* < ta at low Pe, leading to the mixture 
of the diffusion-ballistic-diffusion regimes. 

A rather good fit, not shown here, of the finite density data in the limit 
Pe 3> 1 and for time-delays such that t > t* is achieved by using the single 
dumbbell expression in Eq. (fTOj) 




1 — e 


D 


R 


Pe > 1 , 


(23) 


without the first term (negligible if Pe 3> 1) and upgrading the remaining 
parameters, and 77^, to be density-dependent fitting parameters, as done 


m 


l4 l25| . The quality of this fit was discussed in 


The crossover time-delay between the last ballistic or super-diffusive, and 
the diffusive regimes increases, though rather weakly, with 0, at low Pe, see 
the inclined dashed line in the last panel that is also a guide-to-the-eye. 

In 24 we showed that the qualitative dependence of Da on ksT for the 
single dumbbell case (0 = 0), is maintained under interactions (0 7 ^ 0). For 
Pe -C 1, Da is dominated by thermal fluctuations and increases with ksT. 
Instead, for Pe 1, Da is dominated by the work done by the active forces. 
Da saturates at small values of T for 0 > 0.2. Instead, at high temperatures 
Da seems to retain the linear growth with temperature of the single dumbbell 
at least for the temperatures used in the simulations. 

The 0 dependence of Da at fixed T and for different active forces was 
discussed in 23, 24| where it was shown how the Tokuyama-Oppenheimer 
law of the passive system was simplified under activation. On the one hand, 
we showed that the ratio of diffusion coefficients of the active system at finite 
density and single passive dumbbell 77^(Fact, T, 0)/77 p^ depends on Fact and 
T only thr oug h Pe, as it does for the single dumbbell. The numerical data 
suggested 24 


F^(Fact, T, 0 ) ~ Dli (1 + PeVS) (24) 

with b a non-monotonic function of Pe for Pe = 40 and Pe = 66 , and 
DAiFsu^t, T, 0) ~ Fa (Fact, T, 0) [1 ai(Pe) 0 -f 03 (Pe) 0^] (25) 
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for Pe = 4 and Pe = 20, with oi negative in all cases while 02 changing 
sign from negative at Pe < 20 to positive at Pe > 20 (leading to a growing 
behaviour at large (j) that is not physical). At Pe = 20 the density dependence 
is almost linear as 02 is very close to zero. 




Figure 3: Normalized distribution of the centre of mass horizontal displacements, Aa:, 
in a system with density cj) = 0.1. Four Peclet numbers are used in each panel, Pe = 
40, 20, 4, 2 or, equivalently, T = 0.005, 0.01, 0.05, 0.1, as indicated in the keys. The 
time-delay in each of the four panels is t = t/, 10 t/, ta/2, 10 to (from left to right and 
from top to bottom). ti = 0.1 and the values of to depend on Pe as given in Table [T] The 
probability distribution functions are scaled with ax = The solid curves are 

Gaussian pdfs in normal form (with zero mean and unit variance) and represent the data 
in the first two panels (t = t/, 10 tj) rather accurately. Instead, the data for high Pe in 
the last two panels (t = ta/2, 10 to) are not Gaussian. See the main text for a discussion. 
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4-2. The fluctuations of the centre of mass displacements 

We study the distribution of centre of mass displacements 

1 ^ 

F(Aa:) = — ^(5 (Ax - {Xcraft + tfl - Xcrafltfl))) (26) 

i=l 

or the self-part of the van Hove correlation function. 

In the zero density limit and under no active force the distribution of Ax 
is Gaussian. We will now determine how this limiting form is modihed by 
the active force and the interactions in the various time-delay regimes. 

In Fig. [3] we display the statistics of the horizontal center of mass dis¬ 
placements in a system with a relatively low density, 0 = 0.1. Data for Pe = 
40, 20, 4, 2 are gathered in each of the panels that correspond to different 
time-delays, t = f/, 10 f/, fa/2 and 10 ta- The data are scaled by the hor¬ 
izontal contribution to the mean-square displacement, ax = (Ax^)^/^. The 
solid lines are Gaussian pdfs in their normal form and describe the data for 
short time-delays and all Pe in the hrst two panels very accurately. Instead, 
the data in the last two panels are close to the Gaussian form only for the 
low Pe’s while for high Pe’s the shape of the pdfs is different. In the third 
panel in Fig. |3]the time-delay, t = ta/2, is such that the mean-square dis¬ 
placements are in the second ballistic regime, that in the single dumbbell 
case corresponds to a^ = (Ar^^^) ~ (Fact/q)^ The distribution of Ar is 
peaked around a constant value of Factfa/( 27 ) which projected on one axis 
gives this double shape. The down-pointing arrows inside the plot indicate 
these instants for Pe = 40 and are close to the location of the peaks in the 
hnite density case. 

From the Gaussian £t of the bare data for F(Ax) in the last diffusive 
regime (f = 10 ta, fourth panel in Fig. [3]) we extracted Da and we found 
good agreement with the values of Da obtained in the analysis of the mean- 
square displacement (not shown). 

In the following two hgures. Figs. |4]and|5l we checked the density depen¬ 
dence of the centre of mass displacement distribution. From the hrst to the 
fourth panels t = tj, 10 tj, ta/^,, 10 ta, respectively. In each panel data for 
(j) = 0.01, 0.1, 0.3, 0.5, 0.7 are shown. 

In the cases in which activation is strong. Fig. 01 at short time delays 
(hrst two panels) all systems have Gaussian huctuations. At long-time delays 
(last two panels) the distributions depend on </. In the regime III there is 
a two peak structure at low density, while the central part becomes close to 
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Figure 4: Distribution of center of mass translation displacement in the horizontal direc¬ 
tion, Ax, for Pe = 40 (or T = 0.005), and various densities if) = 0.01, 0.1, 0.3, 0.5, 0.7. 
The four panels correspond to time-delays in the four regimes labelled I, II, III and IV 
(t/, 10 ti, ta/2 and 10 ta, respectively). For I and II the Gaussian in normal form de¬ 
scribes the data very well. In the last two panels the statistics are not Gaussian and there 
is a non-trivial density and time-delay dependence that we discuss in the main text. The 
straight lines in the third panel {t = ta/2) are exponential fits to the tails of the (p = 0.7 
data. 


flat at higher densities, due to stronger interactions, and then develops two 
exponential wings at larger absolute displacements. In the regime IV the 
distributions do not have the double peak structure and are Gaussian close 
to Ax ~ 0 for all densities. Looking carefully at the MSD for Pe = 40 in the 
left panel in Fig. |2l one notices that the curves in this regime are less flat 
at low density while they are closer to being flat at 0 = 0.5. The crossover 
to normal diffusion is, therefore, slower for lower density. The long-lasting 
super-diffusive behaviour seems to be related to the non-Gaussianity of the 
pdfs. The MSD for Pe = 2 in the right panel in Fig. [2] are flatter earlier for 
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cf) < 0.5 and, consistently, the pdfs are Gaussian for all these densities (see 
Fig. [5]). The case 0 = 0.7 is different since the system is far from the last 
diffusive regime for the time-delays used. 

In the cases in which activation is weak. Fig. |5l once again at short 
time delays (hrst two panels) all systems have Gaussian fluctuations. At the 
intermediate time delays (third panel) the distributions depend on 0 and, 
for high densities, see (j) = 0.5, 0.7 in the plot, there is an excess weight on 
large deviations with respect to the Gaussian, that is close to exponential. 
The effect becomes more evident when the sub-diffusive behavior is more 
pronounced, see Fig.|2l At still longer time-delays (fourth panel) the statistics 
becomes Gaussian again for all densities, except for the case at 0 = 0.7 since 
in this case the time-delay is still too short to reach the hnal diffusive regime. 


5. Rotational motion at finite density 


We turn now to the rotational dynamics. In Fig. |6]we display the angular 
MSD normalized by time-delay. The two panels show data for Fkct = 0.1 and 
Pe = 40 and Pe = 2, and the same densities as in Fig. |2J These plots also 
show interesting features: 

- In all cases there is a hrst ballistic regime with a pre-factor that is inde¬ 
pendent of 0 and increases with temperature (The case t tj = rud/y of 
the single dumbbell.) 

- Next, the dynamics slow down and, depending on T and 0, the normalised 
MSD may attain an ever-lasting plateau associated to normal diffusion for 
low 0 at any temperature, or even decrease, suggesting sub-diffusion, at high 
enough 0. 

- At high Pe and sufficiently high density the dynamics accelerate next, with 
a second super-diffusive regime that crosses over to a hnal dihusive regime. 

- In the late normal dihusive regime all curves saturate and the height of the 
plateau yields the diherent Dr coefficients that we discuss below. 

In the phase separated regime the dumbbell clusters rotate (3, [ 23 . It 
is possible that strong huctuations not far from the critical point have an 
important rotational component than enhances/advects rotational dihusion 
giving rise to an observable contribution to displacement. 

The dihusion constant decreases with temperature below a crossover be¬ 
yond which it increases approximately linearly with temperature 2^. The 


Fact-independence of Dr is lost as soon as the interaction between dumbbells 
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Figure 5: Distribution of the center of mass translation displacement in the horizontal 
direction, Ax, for Pe = 2, and various densities (j) = 0.01, 0.1, 0.3, 0.5. The four panels 
correspond to time-delays t = tj, 10 ti, tal2 and 10 ta, respectively. The first and last 
cases are in the regimes labelled 1 and IV for (j) < 0.5. As in Fig. |4] the tails of the 
distribution for (p = 0.7, and t = tal^ and t = IQ ta in this case, are well described by 
exponential fits, represented as straight lines. See the main text for a discussion. 


is switched on. The density dependence is also quite complex and was dis¬ 
cussed in j^. Finally, we showed that Dji/{kBT) depends on Fact and ksT 
only through Pe, 

F^(Fact,T,0) = fcBT/R(Pe,0) (27) 

with /R(Pe,0) = /j?(0,0) = 2 /( 70 -^). At low densities, while the master 
curve decreases with 0 in the whole range considered for Pe < 20, it becomes 
flat at Pe = 20 and it increases with 0 for Pe >20, for cp < 0.5. This would 
suggest: 

/it(Pe, 0) ^ ^ + a(Pe, (p) , (28) 
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Figure 6: The angular MSD for Pe = 40 and Pe = 2 and, with different lines, various 
densities given in the keys to each panel. The dash in the first panel highlights the initial 
ballistic behavior. The vertical arrows indicate the characteristic times t/, t*, ta- Note 
that the vertical scale is different in the two panels. 


with a(Pe, 0) almost linear in 0 and the slope changing sign at Pe ~ 20 for 
small (j). At all Pe there is a cross-over at high enough densities after which 
the rotational diffusion constant decreases with increasing density. One can 
associate this feature to the fact that for sufficiently dense systems rotations 
are inhibited and Dji decreases. 

5.1. The angular displacement distribution 

Figure [7] demonstrates that for low enough density, 0 = 0.1 in this hgure, 
the angular displacement distribution function at all time delays and for all 
Peclet numbers can be put into a normal Gaussian form after normalisation 
by ag = 

In Figs. [H] and [9] we analyse the density dependence of the normalised 
probability distribution of angular displacements. In Fig. |8] a high Peclet 
number is used, Pe = 40, while in Fig. [9] the Peclet number is low, Pe = 2. 
In the hrst case, the four time-regimes I, II, III, IV, in the behaviour of the 
center of mass translational mean-square displacement, are well separated. 
In the second case, they are not. We see that for 0 = 0.01, 0.1, 0.3 the data 
in both hgures collapse onto the normal Gaussian. The higher density data, 
0 = 0.5, deviate from this master curve when the time delay is chosen to 
be f = 10 ti for Pe = 40, and t = G/2 = 25 for Pe = 2. The data points 
for 0 = 0.7 in the last two panels in the two hgures have a peak at A9 = 0 
corresponding to dumbbells that do not rotate between the two times and 
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Figure 7: Normalized distribution of angular displacements, A9, in a low-density system, 
(p = 0.1. Four Peclet numbers are used in each panel, Pe = 40, 20, 4, 2, as indicated in 
the keys. The time-delays in the four panels are t = t/, 10 ti, ta/2, 10 ta- The values 
of ta depend on Pe and are given in Table [T] The probability distribution functions are 
scaled with ag = and the solid lines are a Gaussian in normal form, with zero 

mean and unit variance. 


two exponential wings ty pica l of heterogeneous systems. Note that these 
pdfs resemble the ones in ^ though for translational displacements in this 


case. 


6. Conclusions 

We considered the diffusion properties of a system of active dumbbells 
with repulsive interaction. We focused on the Peclet number and density 
regime in which the global system is homogeneous and we studied the fluc¬ 
tuations of the centre of mass translation and angular displacements. 

We hrst summarised the translational and rotational MSD of the system. 
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Figure 8: Normalized distribution of angular displacements A9 for Pe = 40. Data for 
various densities (p = 0.01, 0.1, 0.3, 0.5, 0.7 are collected in each panel. In each of the 
four panels the time-delay {t = ti, 10 t/, ta/2, 10 ta) is chosen to lie in one of the four 
relevant time-scales called I, II, III, IV. t* = 10 <C ta = 1000, see Table [U for this high 
value of Pe. The solid curve is a Gaussian in normal form. In the second, third and fourth 
panels deviations from the Gaussian are visible for the high density data, (p = 0.7. In the 
last two panels the straight lines are exponential fits to the tails. 


In the single particle limit, the translational MSD has a very rich time-delay 
dependence, with four distinct time regimes (ballistic, diffusive, ballistic and 
diffusive) separated by three characteristic times (the shortest inertial, ti, 
the diffusive, ta, and an additional one, t* oc ta/Pe^, that lies in between the 
other two for large Pe). This rich structure survives under hnite densities 
with modihed parameters. The diffusion constant in the last diffusive regime 
has a non-monotonic dependence on temperature, as for the single dumbbell 
case, and it decreases with increasing self-propelled particle density at all 
temperatures. Moreover, it depends on temperature and active force only 
through the Peclet number at all densities explored. In general, the Pe 


18 









Figure 9: Normalized distribution of angular displacements A6 for Pe = 2. Data for 
various densities (j) = 0.01, 0.1, 0.3, 0.5, 0.7 are collected in each panel. The time-delay 
in the first {t = ti) and last {t = 10 fa) panels lie in the inertial ballistic and last diffusive 
regimes I and IV, respectively, for (j) < 0.5. The time-delays in the second and third panels 
(t = 10 ti and t = ta/^, respectively) are in the intermediate crossover between these 
two. t* = 4000 '> ta = 50, see the Table [H for this low Pe number. The solid curves 
are Gaussians in normal form. The deviations from the Gaussian in the intermediate 
regime are pronounced for lower densities than for Pe = 40. The solid straight lines are 
exponential fits to the tails in the (j) = 0.7 data. 


dependence is non-nionotonic. 

The behavior of the rotational MSD, that is rather simple for a single 
dumbbell with just one crossover between ballistic and diffusive behaviour, 
reflects the behavior ot the translational MSD at hnite densities, where inter¬ 
mediate regimes also appear. The late epochs rotational diffusion constant 
increases with temperature (though not linearly) at all densities and active 
forces simulated and depends on temperature and activity only through Pe. 
At low densities, its dependence on density changes from decreasing at low 
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Pe to increasing at high Pe. 

We have then evaluated the distribution functions for the centre of mass 
and angular displacements for time-delays corresponding to the various dy¬ 
namical regimes. For time-delays shorter or of the order of the inertial time 
tj the distributions were always found to be very close to a Gaussian, for 
all the (j) and Pe considered. At large Pe, always in the homogeneous phase 
but not far from the critical value (Pe/Pec ~ 0.62) j^, for time-delays 
corresponding to the super-diffusive or second ballistic regime of the single 
particle, while the angular distributions remain Gaussian except at the high¬ 
est density considered 0 = 0.7, the translational displacement distributions 
show a two-peak character at low density becoming more rounded when the 
density increases. The position of the peaks corresponds to the analytic esti¬ 
mate given in Sec. 13.11 The distributions get closer to a Gaussian in the last 
regime. The effects on the final diffusive regime due to the cross-over with 
the preceding super-diffusive regime are more pronounced than at small Pe 
even at times of the order of 10 ta- Accordingly, non-Gaussian effects can be 
seen in the translational distributions at these times. The case at 0 = 0.7 is 
different since, for the time-delay considered, the system is still far from the 
hnal diffusive regime. At small Pe, when the effects of interactions are less 
relevant (Pe/Pec ~ 0.03), the second regime becomes sub-diffusive when the 
density increases. In correspondence with this, deviations from the Gaussian 
are found in both translational and angular displacement distributions. Fi¬ 
nally, both distributions have a Gaussian character in correspondence with 
the last diffusive regime, except at </ = 0.7. In this very high density limit we 
hnd exponential tails in the translational and rotational distributions, and a 
central peak at vanishing angular displacement. 

After this work we plan to analyse the motion of tracers in contact with 
this active sample and, especially, to analyse the existence of a parameter to 
be interpreted as an effective temperature from the mobility and diffusive 
properties of the sample and the tracers, as done in 28, 301 for another 

active matter model and in 3l| for a sample of active Janus particles. 
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